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Abstract. 

We suggest to use the observationally measured and theoretically justified correlation between size and rotational velocity of 
galactic discs as a viable method to select a set of high redshift standard rods which may be used to explore the dark energy 
content of the universe via the classical angular-diameter test. Here we explore a new strategy for an optimal implementation of 
this test. We propose to use the rotation speed of high redshift galaxies as a standard size indicator and show how high resolution 
multi-object spectroscopy and ACS/HSF high quality spatial images, may be combined to measure the amplitude of the dark 
energy density parameter Qg, or to constrain the cosmic equation of state parameter for a smooth dark energy component 
(w = pip, -1 < w < -1/3). Nearly 1300 standard rods with high velocity rotation in the bin V = 200 ± 20km/s are expected 
in a field of 1 sq. degree and over the redshift baseline < z < 1.4. Fhis sample is sufficient to constrain the cosmic equation 
of state parameter w at a level of 20% (without priors in the [Cl,„, Q.q] plane) even when the [OID/13727A linewidth-diameter 
relationship is calibrated with a scatter of ~ 40%. We evaluate how systematics may affect the proposed tests, and find that a 
linear standard rod evolution, causing galaxy dimensions to be up to 30% smaller at z = 1.5, can be uniquely diagnosed, and will 
minimally bias the confidence level contours in the [CIq, w] plane. Finally, we show how to derive, without a priori knowing 
the specific functional form of disc evolution, a cosmology-evolution diagram with which it is possible to establish a mapping 
between different cosmological models and the amount of galaxy disc/luminosity evolution expected at a given redshift. 

Key words, cosmology: observations — cosmology:fheory — cosmologyxosmological parameters — galaxies: distances and 
redshifts — galaxies: fundamental parameters 



1. Introduction 

Several and remarkable progresses in the understanding of the 
dynamical status of the universe, encourage us to believe that, 
after roaming from paradigm to paradigm, we are finally con- 
verging towards a well-founded, internally consistent standard 
model of the universe. 

The picture emerging from independent observations and 
analysis is sufficiently coherent t o be r eferred to as the con- 
cordance model (e.g. Tegmark, 2006). Within this frame- 



work, the universe is fiat (£2# = -O.OOS^qqjqj) composed 



of ~ 1/5 cold dark matter (Q, 
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dark energy (£2a = 



0.761 tj5 nil), with large negative pres- 



sure (w = -0.941^QQ}g), and with a very low baryon content 
(Qb - 0.0416jqqi|). Mounting and compelling evidence for 
accelerated expansion of the universe, driven by a dark en- 
ergy component, presently relies on our comprehension of the 
mec hanisms with wh i ch Sup e rnovae la (SNIa) emit radiation 



(see Perlmutter et al. (1999); IRiess et al.1 (120011) ) and of the 



physical processes that p roduced t e mpera t ure fluctuations in 
the pr i meval plasma (se e Lee et al.l J200lb: de Bernardis et al 



(2002T): lHalverson et al.1 d2002l) :ISpergel et al. (2006).) 



Even if the ambitious task of determining geometry and 
evolution of the universe as a whole, which commenced in the 
1930s, now-day shows that the relativistic Friedman-Lemaitre 
model passes impressively demanding checks, we are faced 
with the challenge of developing and adding new lines of 
evidence supporting (or falsifying) the concordance model. 
Moreover, even if we parameterize our ignorance about dark 
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energy describing its nature only via a simple equation of state 
w = pip, we only have loose constraints on the precise value 
of the w parameter or on its functional behavior. 

In this spirit we focus this analysis on possible complemen- 
tary approaches to determining fundamental cosmological pa- 
rameters, specifically on geometrical tests. 

A whole arsenal of classical geometrical methods has been 
developed to measure global properties of the universe. The 
central feature of all these tests is the attempt to directly probe 
the various operative definitions of relativistic distances by 
means of scaling relationships in which an observable is ex- 
pressed as a function of redshift (z) and of the fraction of criti- 
cal density contributed by all forms of matter and energy (Q). 

The most remarkable among these classical methods are 
the Hubble diagram (or magnitude-redshift relation m - 
m(M, z, Q)), the Angular diameter test (or angle-redshift re- 
lation 9 - 9{L, z, fi)), the Hubble test (or count-redshift rela- 
tion N=N(n,z,Q)) or the Alcock-Paczinsky test (or deformation- 
redshift relation Az/zA# = k — k(z, fi)). The common key idea 
is to constrain cosmological parameters by measuring, at var- 
ious cosmic epochs, the scaling of the apparent values m, 0, 
N, k of some reference standard in luminosity (M), size (L), 
density (n) or sphericity and compare them to corresponding 
model predictions. 

The observational viability of these theoretical strategies 
has been remarkably proved b y the Supernova Cosmology 



Project ( Perlmutte r et al. , 1999 ) and the High-z Supernova 



Team dRiess et all 2001) in the case of the Hubble diagram. 
With a parallel strategv. lNewman et alJd2002l) recently showed 
that a variant of the Hubble test (N(z) test) can be in prin- 
ciple applied to distant opti cal clusters selected in deep red- 



shift survey such as VVDS dLe Fevre et all 120051) and DEEP2 



dDavis et al. , 2000), in order to measure the cosmic equation- 
of-state parameter w. 

Unfortunately, the conceptually simple pure geometrical 
tests of world models, devised to anchor relativistic cosmology 
to an observational basis, have so far proved to be difficult to 
implement. This is because the most effective way to constrain 
the evolution of the cosmological metric consists in probing 
deep regions of the universe with a primordial class of cosmo- 
logical objects. Besides the complex instrumental technology 
this kind of experiments requires, it becomes difficult at high 
redshift to disentangle the effects of object evolution from the 
signature of geometric evolution. 

Since geometrical tests are by definition independent from 
predictions of theoretical models or simulations, as well 
as from assumptions about content, quality and distribution 
of matter in the universe (mass fluctuations statistics (e.g. 
Haiman, Mohr & Holder 2001, Newman et al. 2001), galactic 
biasing (e.g. Marinoni et al. 1998, Lahav et al. 2001, Marinoni 
et al. 2006), Halo occupation models (e.g. Berlind et al. 2001, 
Marinoni & Hudson 2002, van den Bosch et al. 2006)) it is of 
paramount importance to try to devise an observational way to 
implement them. The technical maturity of the new generation 
of large telescopes, multi object spectrographs, large imaging 
detectors and space based astronomical observatories will al- 
low these tests to be mor e effectively applied in the near future 
(IHuterer & Turnen l2000l) . In this paper, we describe a method 



to select a class of homologous galaxies that are at the same 
time standard in luminosity and size, that can be in principle ap- 
plied to data coming from the zCOSMOS spectro-photometric 
survey (Lilly et al. 2006, of the deep universe. 

An observable relationship exists between the speed of 
rotation V of a spiral galaxy and its metri c radial dimen- 
sion D as well as its tota l luminosity L dTullv & Fishei , 
1977 ; Bottinelli et al. , 1980). From a theoretical perspective, 
this set of scaling relations are expected and explicitly pre- 
dicted in the context of CDM models of galaxy formation 
dMo, Mao & Whitel Il998h . The Tully-Fisher relations for di- 
ameter and luminosity have been extensively used in the local 
universe to determine the distances to galaxies and the value of 
the Hubble constant. We here suggest that they may be used in 
a cosmological context to select in a physically justified way, 
high redshift standard rods since galaxies having the same rota- 
tional speed will statistically have the same narrow distribution 
in physical sizes. 

The picture gets complicated by the fact that the standard 
model of the universe implies some sort of evolution in its con- 
stituents. In a non static, expanding universe, where the scale 
factor changes with time, we expect various galaxy properties, 
such as galaxy metric dimensions, to be an explicit function 
of redshift. In principle, one may break this circular argument 
between model and evolution with two strategies: either by un- 
derstanding the effects of different standard rod evolutionary 
patterns on cosmological parameters, or by looking for cosmo- 
logical predictions that are independent from the specific form 
of the disc evolution function. 

In this paper, following the first approach, we study 
how different disc evolution functions may bias the angular- 
diameter test. We simulate the diameter-redshift experiment us- 
ing the amount of data and the realistic errors expected in the 
context of the zCOSMOS survey. We then evaluate how differ- 
ent disc evolution functions may be unambiguously recognized 
from the data and to what extent they affect the estimated val- 
ues of the various cosmological parameters. We also explore 
the second approach and show how cosmological information 
may be extracted, without any knowledge about the particular 
functional form of the standard rod evolution, only by requiring 
as a prior an estimate of the upper limit value for the relative 
disc evolution at some reference redshift. 

This paper is set out as follow: in §2 we review the the- 
oretical basis of the angular diameter test. In §3 we describe 
the proposed strategy to select high redshift standard rods. In 
§4 we digress on how implement in practice the '6 — z' test 
with zCOSMOS data, and in §5 we present the zCOSMOS ex- 
pected statistical constraints on cosmological parameters. In §6 
we discuss different possible approaches with which to address 
the problem of standard rod evolution. Conclusions are drawn 
in §7. 

2. The Angular Diameter Test 

We investigate the possibility of probing the cosmological met- 
ric using the redshift dependence of the apparent angular diam- 
eter of a cosmic standard rod. What gives this test special ap- 
peal is the possibility of detecting the "cosmological lensing" 
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effect, which causes incremental magnification of the apparent 
diameter of a fixed reference length. 



n„=o.3 n =o.7 w=-i.o 



Let's consider the transverse comoving distance (see Hogg, 

1999) 



r(z, n m , Q Q , w) = C -==S k { VT^~i f E(xT l dx) 

H V Qt Jo ' 



where 



E(x) = [Q,„(l + x)" + Q e (l + x) 3+3w + Q*(l + xY] 



y2 n l/2 



(1) 



(2) 



and where S -i(y) = sinh(y), S i(y) = sin(y), So(y) = y while 
— 1 — £l m — £1q. 

An object with linear dimension D at a redshift z has thus 
an observed angular diameter 9 



0(z,p) 



D 



r(z,p) 



(1+z) 



(3) 



which depends on the general set of cosmological param- 
eters p = [Q. m , Q.q, w] via the relativistic definition of angular 
distance, d\ = r(z,p)/(l +z). 

This test may be implemented without requiring the knowl- 
edge of the present expansion rate of the universe (the depen- 
dence from the Hubble constant cancels out in eq.0. At vari- 
ance, although characterized by a smooth and diffuse nature, 
dark energy significantly affect the dynamic of the universe. 
From eq. Q] it is clear that the angular-diameter test depends 
on the dark energy component via the expansion rate evolution 
E(z). The more negative w, the more accelerated the expansion 
is and the smaller a fixed standard rod will appear to an ob- 
server. 

The efficiency of different cosmological observables in 
probing the nature of space-time ultimately depends upon their 
sensitivity to the cosmological parameters Q,„, Q.q, w. The rel- 
ative sensitivity of empirical cosmological tests based on the 
scaling of the angular diameter distance (aU) and of the volume 
element (dV/dz = (c / Ha)(r 2 / E(z)) is derived in Fig.Q] where 
we assume that Poissonian errors are constant in time and no 
redshift-dependen t syste matics perturb the measurements (e.g. 
Huterer & Turner! §000)). Since the luminosity distance (i.e. 
the distance inferred from measurements of the apparent mag- 
nitude of an object of known absolute luminosity) is defined 
as di = (1 + z) 2 dA, we note that the angular diameter test has 
the same cosmological discriminatory power as the Hubble di- 
agram. The upper panel of Fig. Q] shows that the sensitivity of 
both d& and dV/dz to the mean mass density parameter, Q m , 
increases monotonically as a function of redshift. This means 
that the deeper the region of the universe surveyed, the more 
constrained the inferred value of Q m is. 

Conversely, the sensitivity of both empirical tests to a 
change in the constant value of w peaks at redshift around unity, 
and levels off at redshifts greater than ~ 5. The reason for this is 
that the dark energy density pq, which substantially contributes 
to the present-day value of the expansion rate was negligible in 
the early universe (pq/pm 00 (1 + z) 3w , see eq.|2|. 

The fact that we are living in a special epoch, when two 
or more terms in the expansion rate equation make compara- 
ble contributions to the present value of E(z), can be appreci- 
ated in the central panel of Fig. [T] Because each of the terms 
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Fig. 1. The relative sensitivity of the angular diameter distance 
(aU) and volume element (dV/dz) to a change in the values of 
Q. m , Qq and w. The partial derivatives are computed with re- 
spect to the position (Q„, = 0.3, Q.q = 0.7, w = -1) in the 
parameter space. 

in eq. [2] varies with cosmic time in a different way, there is a 
redshift window where the search for Q G is less efficient (i.e. 
2 < z < 4). Therefore one can maximize the cosmological 
information which can be extracted from the classical tests of 
cosmology, and specifically from the angular diameter test, by 
devising observational programs probing a large field of view 
in the redshift range < z < 2. 

3. The Standard Rod 

A variety of standard rod candidates have been explored in pre- 
vious attempts of implementing the angular diameter-redshif t 
test: gal axies dSandagel. Il972l lDiorgovski & Spinradi 1 198 1 ), 



clusters dHicksod| 1977; Bruzual & Sp inradi 1 1 9781; iPenl ll 9971) . 
halo clustering ( Coorav et all 200 lb . Those methods failed to 
yield conclusive evidence because the available redshifts were 
few and local, and the quality of the imaging data used in the 
estimate of sizes was poor. 

Good quality size measurements for high r edshift object s 
have be c ome a vailable for radio sources (e.g. Miley (1971); 
i Kapahil d 19751). and r ecentl y several authors dKellermannl 



1993; Wilki nson et all 119981) have reported a redshift depen- 



dence of radio source angular sizes at 0.5 < z < 3, which is not 
easily reconciled with other recent measurements of the cosmo- 
logical parameters (but see Daly & Djorgovski 2004 for results 
more consistent with the concordance model.) 

The radio source results may be affected by a variety of 
selection and evolutionary effects, the lack of a robust defini- 
tion of size, and by difficulties in assembling a large, homo- 
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geneous sample of radio observati ons dBuchalter et all 1 1998 



Gurvits, Kellermann & Frey, 1999). 



CALIBRATOR 



A common thread of weakness in all these studies is that 
there are no clear criteria by which galaxies, clusters, extended 
radio lobes or compact radio jets associated with quasars and 
AGNs should be considered universal standard rods. Moreover, 
lacking any local calibration for the metric size of the standard 
rod, the standard rod dimension (parameter D in eq. [3]) is often 
considered as a free fitting parameter. Since the inferred cos- 
mological para meters heavily depend o n the assumed value for 
the object size (Lima & Alcanizll2002l) . an a priori independent 
statistical study of the standard rod distribution properties is an 
imperative prerequisite. 

We thus propose to use information on the kinematics of 
galaxies, as encoded in their optical spectrum, a) to identify 
in an objective and empirically justified way a class of ob- 
jects behaving as standard rods, and b) to measure the abso- 
lute value of the standard rod length. The basic idea consists in 
using the velocity-diameter relationship for d isc galaxies (e.g. 



Tully & Fisherl 11977b ISaintonge et all |2007) as a cosmologi 



cal metric probe. In Fig. [2] we plot t he local relationship w e 



have derived in paper II of this series ((Sai ntonge et all 120071) ') 
between half-light diameters and rotational velocities inferred 
using the H„/l6563Aline. The sample used to calibrate the 
diameter-velocity relationship is the SFI++ sample described 
by Springob et al 2007. Also shown is the the amplitude of the 
scatter in the zero-point calibration of the standard rods. 

In the local universe, rotation velocities can be estimated 
from either 21 cm HI spectra or from the H„/16563A optical 
emission lines. However, the H a line is quickly redshifted into 
the near-infrared and cannot be used in ground-based optical 
galaxy redshift surveys at z > 0.4, while HI is not detectable 
much past z = 0.15. Only [OII]/l3727A line widths can be 
successfully used in optical surveys to infer the length of a 
standard rod dimension D at z ~ 1 . Clearly one could obtain 
rotational velocity information for high redshift objets by ob- 
serving the Hq. line with near-IR spectrographs. However it is 
much easier to get large samples of kinematic measurements 
using Oil and multi-slit devices, rather than get sparser sam- 
ples using a single-object, near-IR spectrograph. 

A detailed study of the kinematical information encoded 
in the [Oil] line are presented in Paper II. In that paper we 
have explored the degree of correlation of optical H a and [Oil] 
rotational velocity indicators, i.e. how well the rotation veloc- 
ities extracted from these different lines compare. Moreover, 
we have derived a local diameter-velocity relationship, and we 
have investigated the amplitude of the scatter in the zero-point 
calibration of the standard rods. 

We finally note that the present-day expansion rate sets the 
overall size and time scales for most other observables in cos- 
mology. Thus, if we hope to seriously constrain other cosmo- 
logical parameters it is of vital importance either to pin down 
its value or to devise //((-independent cosmological tests. Note 
that, given the calibration of the diameter-linewidth relation in 
the form Hq = f(V), the 0-expression in equation [3] is effec- 
tively independent of the value of the Hubble constant. 
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Fig. 2. The diameter vs. velocity relationship calibrated in pa- 
per II is plotted using a black line. Shadowed regions represent 
the 1 and 2cr uncertainties in the calibrated relationship. D rep- 
resents the corrected (face on) half light diameter while the ve- 
locity has been measured using the H a line. Dotted lines repre- 
sent the upper and lower relative uncertainty [cr D /D],„, = 0.15 
in diameters. The conservative relative dispersion in the rela- 
tionship assumed in this study ([cr D /D],„, = 0.4) is also pre- 
sented using dashed lines. With this conservative choice we 
take into account that sizes and velocity are measured in the 
high redshift universe with greater uncertainties. Within the in- 
terval centered at V = 200 + 20 km s~' (mean physical galaxy 
dimension ~ 10 + 1.5/T 1 kpc which roughly corresponds to an 
isophotal diameter D25 ~ 20 + 3h~ l kpc), we will select galax- 
ies with total absolute I magnitude Mi - 5 logh = -22.4 (see 
paper II). The velocity selected standard rods at z ~ 1 are thus 
well within the visibility window of deep gala xy redshift sur- 
yeys such as the V VPS dLe Fevre et all 120051) or zCOSMOS 
dLillv et all 120061) . For exemple the VVDS is flux limited at 

/ < 24 and selects objects brighter then Mj 20 + 5 log h at 

z = l. 



4. Optimal Test Strategies 

In this section we outline the optimal observational strategy re- 
quired in order to perform the proposed test. With the proposed 
selection technique, the photometric standard rod D is spec- 
troscopically selected and the sample is therefore free from 
luminosity-size selection effects, that is from the well known 
tendency to select brighter and bigger objects at higher red- 
shifts (Malmquist bias) in flux-limited samples. However it is 
crucial that a large sample of spectra be collected, in order 
to obtain V^V gain over the intrinsic scatter in the calibrated 
V(OII)-diameter relationship. 
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4.1. Galaxy sizes measurement 

Since galaxies do not have sharp edges, their angular diameter 
is usually defined in terms of isophotal magnitudes. However 
since surface brightness is not constant with distance, the suc- 
cess in performing the experiment revolves a round the use o f 
metric rather than isophotal galaxy diameters dSandagelll995l) . 

A suitable way to measure the photometric parameter 6, 
without making any a-priori assumption about cosmological 
models, consists in adopting as the standard scale length esti- 
mator either the half-light radius of the galaxy or the ^-function 
of Petrosian (1976). The Petrosian radius is implicitly defined 



(4) 



i.e. as the radius 6 at which the surface brightness averaged 
inside is a predefined factor tj larger than the local surface 
brightness at 8 itself. 

Both these size indicators are independent of K-correction, 
dust absorption, luminosity evolution (provided the evolution- 
ary change of surface brightness is independent of radius), 
wave-ban d used (if there is no color gradient) and source light 
profile dDiorgovski & Spirirad , 1981 ). 



4.2. Standard rods optimal selection 

The choice of the objects for which the velocity parameter V 
and the metric size is to be measured is a compromise between 
the observational need of detecting high signal-to-noise spec- 
tral and photometric features (i.e. selecting high luminosity and 
large objects) and the requirement of sampling the velocity 
distribution function (n(V)dV ~ V~ 4 for galaxy-scale halos) 
within an interval where the rotator density is substantial. 

Given the estimated source of errors (see next section), and 
the requirement of determining both Qg and w with a preci- 
sion of 20%, we find, guided by semi analyti cal models predict- 
ing the redshift distributi o n of r otators (i.e. Naravan & Whitel 
19881 iNewman & Davisl l2000h . that an optimal choice are 
V = 200 km s _1 rotators. 

In particular, as shown in Paper II, the I band character- 
istic absolute magnitude of the V = 200 km s _1 objects is 
Mj ~ -22.4 + 5 log hjo i.e. well above the visibility thresh- 
old of flux-limited surveys such as zCOSMOS or VVDS (as an 
example, for I a t, = 24 the limiting magnitude of the VVDS, one 
obtains that M/ < -21+5 \ogh-n) at z = 1.5). Therefore, with 
this velocity choice, the selected standards do not suffer from 
any Malmquist bias (i.e. any effect which favors the system- 
atic selection of the brighter tails of a luminosity distribution at 
progressively higher redshifts). 



4.3. The zCOSMOS Potential 

In §2 we have emphasized that an optimal strategy to study the 
expansion history of the universe consists in probing, with the 
angular diameter test, the < z < 2 interval. However, in order 
to collect a large and minimally biased sample of standard rods 



over such a wide redshift baseline, we need the joint availability 
of high quality images and high resolution multi-object spectra. 

Space images allow a better determination of galaxy struc- 
tural parameters (sizes, luminosity, surface brightness, inclina- 
tion etc) at high redshift. In particular, the ACS camera of the 
Hubble Space Telescope can survey large sky regions with the 
key advantages of a space-based experiment: diffraction limited 
images no seeing blurring, and very deep photometry. Ground 
multi-object spectrographs operating in high resolution mode 
allow a better characterisation of the gravitational potential- 
well of galaxies, facilitating the fast acquisition of large sam- 
ples of standard rods. For example, in the spectroscopic res- 
olution mode R=2500 (1 "slits), the VIMOS spectrograph (Le 
Fevre et al. 2003) allows to resolve the internal kinematics of 
galaxies via their rotation curve or line-width. Note that the 
VIMOS slitlets can be tilted and aligned along the major axis 
of the galaxies in order to remove a source of potentially sig- 
nificant error in the estimate of the rotation velocities. 

These observational requirements are mandatory for a suc- 
cesfull implementation of the proposed cosmological probe. 
The practical feasibility of the strategy is graphically illus- 
trated in Fig [3] In this figure we show a high resolution spectra 
of a galaxy at z=0.5016 obtained with a total exposure of 90 
min with VIMOS. The ground and space (ACS) images of the 
galaxy are also shown for comparison. 

Interestingly, a large sample of rotators can be quickly as- 
sembled by the currently underway zCOSMOS deep redshift 
survey, which uses the VIMOS multi-object spectrograph at the 
VLT to target galaxies with ACS photometry in the 2 sq. degree 
COSMOS field (Scoville et al. 2003)). In principle, one can 
measure in high resolution modality R=2500 the line-widths 
of [OII]/l3727 A up to z ~ 1 .4 by re-targeting objects for which 
the redshift has already been measured in the low-resolution 
(R=600) zCOSMOS survey. Since the spectral interval covered 
in high resolution mode is limited to ~2000A, this allows us to 
determine the optimal telescope and slitlets position angles in 
order to maximize the number of spectra whose Oil emission 
lines fall onto the CCDs. With this fast follow-up strategy one 
could be able to target about 200 rotators per pointing with an 
exposure time of 4 hours down to I a b = 24. 

In the following we will use the zCOSMOS sample as a test 
case to assess the performances of the proposed cosmological 
test. 

5. Constraints on Cosmological Parameters 

In this section we present in detail some merits and advan- 
tages of the proposed approach for constraining the value of 
the fundamental set of cosmological parameters. We evalu- 
ate the potential of the test, within the zCOSMOS operational 
specifications, in placing constraints not only on the simplest 
models of the universe, which include only matter and a cos- 
mol ogical constant, but also on so-called "quinte ssence" mod- 
els dTurner & Whitel 1 19971; Newman et all 120021) . For the pur- 
poses of this study, we assume in this section w to be constant 
in time up to the redshift investigated z ~ 1 .4. 

We assume that a V(OII)-diameter relation can be locally 
calibrated, and that the diameter of V = 200 km s _1 rota- 
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Z=0.5016 
F775W=19.5 
1=39 deg 
PA=43 deg 





Fig.3. Upper: Public release image of the galaxy a = 53.1874858, 5 = -27.910975 at redshift 0.5016, taken with the filter 
F775 W (nearly I band) of the ACS camera by the GOODs program. For comparison the same galaxy as imaged in the EIS survey 
with the WFI camera at the ESO 2.2mt telescope at La Silla. Botton left: raw spectrum of the galaxy taken by VIMOS at the 
VLT-UT3 telescope with a total exposure time of 90 min and a spectral resolution R=2500. Slitlets have been tilted according to 
the major axis orientation (Position Angle = 43°). Botton right: final processed spectrum showing the rotation curve as traced 
by the Hb (^4859A) line. 



tors may be inferred with a worst(/optimal) case relative error 
[o- D ID] in , = 40%(/15%) (see §3 and Fig |2]) The main con- 
tributions to this error figure are uncertainties in measuring 
linewidths, galaxy inclinations and the intrinsic scatter of the 
empirical relation itself. We then linearly combine this intrin- 
sic scatter with the uncertainties with which the Petrosian or 
half-light radius may be determined from the ACS photometry. 
The error on the half light radii is about erg = +0.04" almost in- 
dependent of galaxy sizes in the range 0. 1-1" (S. Gwyn, private 
communication). 

We consider as the observable of the experiment the log- 
arithm of the angle subtended by a standard rod. We use the 
logarithms of the angles rather than angles themselves because 
we assume that the object magnitudes, rather than diameters, 
are normally distributed around some mean value. Moreover, 
in this way the galactic diameter becomes an additive parame- 
ter, whose fitted value (when a z=0 calibration is not available) 
does not distort the cosmological shape of the 9(z) function. 



The observed values log 6 are randomly simulated around 
the theoretical value logff' (cfr. eq.0 using the standard devi- 
ation given by 



\ D Jim \ Jobs 



1/2 



(5) 



For the purposes of this section, 8 r is computed assuming 
a flat, vacumm dominated cosmology with parameters Q,„ 
0.3. Q.q = 0.7, and w = — 1 as reference model. The confidence 
with which these parameters are constrained by noisy data is 
evaluated using the x 1 statistic 



X' 



z 



[Iog0,-log0f(z,/;)] 2 



(6) 



where ff h is given by eq.[3] 

We also derive the expected redshift distribution of galax- 
ies having circular velocity in the 180 < V < 220 km s _I 



200 



150 



z 100 



50 



0.5 



Marinoni et al.: Geometrical tests of cosmological models. I, 

Simulation of the experiment 
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Fig. 4. Left: redshift evolution of the diffe rential comoving numb er density of halos with a circular velocity of 200 km s -1 
computed according to the prescriptions of Newman & Davisl ( 2000l) in the case of a flat cosmological model having Q m = 0.3 
today and a cosmological constant. Only a fraction f=0.2 of the total predicted abundance of halos (i.e. ~ 1300 objects per square 
degree) is conservatively supposed to give line-width information useful for the angular-diameter test. Right: simulation of the 
predicted scatter expected to affect the angular diameter-redshift diagram should in principle achieve with the angular diameter 
test. The simulation is performed assuming the sample is composed by ~ 1300 rotators with V = 200 km s _1 and that a flat 
model with parameters [Q. m = 0.3, Q.q = 0.7, w = -1] is the true underlying cosmological framework. The circular velocity 
ha s been converted into a n estimate of the galaxy diameter (D v = 20 kpc) by using the velocity-diameter template calibrated 
by iBottinelli et al. (1980). The worst-case scenario ([cr D /£)],„, = 40%) is presented. The solid line visualizes the underlying 
input cosmological model 6 ACDM (z), while triangles are drawn from the expected Poissonian fluctuations. The dot-dashed line 
represents the expected scaling of the angular diameter in our best recovered cosmological solution. The dashed and the dotted 
lines represent the angular scaling in a Einstein-de Sitter and Euclidean (non-expanding type cosmology with zero curvature) 
geometry respectively. 



range, in various cosmological scenari os within the framework 
of the Press & Schechte r formalism (INarayan & WhiteL 1 1988 ; 
Newman & Davis , 2000). We note that due to the correlation 
between circular velocity and luminosity, these galaxies could 
be observed to the maximum depth (z ~ 1 .4) out to which the 
Oil is within the visibility window of VIMOS. We take into ac- 
count the uncertainties in the semi-analytic predictions and our 
ignorance about the fraction of discs to be observed that will 
have a spectroscopically resolved [OII]/13727A line, by mul- 
tiplying the calculated halo density by a "conservative" factor 
f=0.2. Using the VIMOS R=2500 resolution mode data, we 
thus expect to be able to implement our angular-diameter re- 
search program using nearly 1300 standard rods per square de- 
gree, which is what has been simulated (see Fig. [4]) 



Since we do not make a-priori assumptions about any pa- 
rameter, and in particular we do not assume a flat cosmology, 
results should be distributed as a xl with v — 3 degrees of 
freedom, which can be directly translated into statistical confi- 
dence contours, as presented in Fig. [5] This figure shows that 
even without assuming a flat cosmology as a prior and consid- 
ering a diameter-linewidth relationship with a 40% scatter, by 
targetting ~ 1300 rotators we can directly infer the presence 
of a dark energy component with a confidence level better than 
3<x. At the same time, its equation of state can be constrained 
to better than 20% (~ 10% if the the diameter-linewidth rela- 
tionship is calibrated with a ~ 15% relative precision). 
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Fig. 5. Predicted 1,2, and 3cr confidence level contours for application of the angular-diameter test. The likelihood contours 
have been derived by adopting a A-cosmology [0,„,0g,w] = [0.3,0.7,-1] as the fiducial model we want to recover, and by 
conservatively assuming one can obtain useful line-widths information for a sample of galaxies having the redshift distribution 
and the Poissonian diameter fluctuations simulated in Fig. [4] Confidence contours are projected onto various 2D-planes of the 
[f2,„, Oq, w] parameter space, and the jointly best fitted value along the projection axis, together with the statistical significance 
of the fit, are reported in the insets. Note the strong complementarity of the confidence region orientation which is orthogonal to 
the degeneration axis of the CMB measurements. Top:) constraints derived assuming that one might survey only 1 square degree 
of sky and that the V(OII)-diameter relation is locally calibrated with a 40% of relative scatter in diameter. Center:) as before but 
assuming a scatter of 15% in diameters. Bottom:) confidence contours for a survey of 16 deg 2 (which corresponds to the full area 
surveyed by VIMOS-VLT Deep Survey) assuming a template V(OII)-diameter relationship with a scatter of 30% in diameters. 



6. Standard Rod Evolution 

The previous analysis shows that the angular diameter test, 
when performed using fast high resolution follow-up of zCOS- 
MOS spectroscopic targets may be used as a promising addi- 



tional tool to explore the cosmological parameter space and di- 
rectly measure a dark energy component. However, the impact 
of any standard rod evolution on these results needs to be care- 
fully examined. 
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Fig. 6. Best fitting cosmological parameters inferred by applying the angular diameter test to data affected by evolution. The 
output (biased) estimates of D.,„ and Oa are plotted as a function of the relative diameter evolution for the following evolutionary 
models: linear (solid line), square-root (dotted line) and quadratic evolution (dashed line). The biasing pattern is evaluated for 
three different fiducial cosmologies: a flat, A-dominated cosmology (£2a = 0.7, left), a low-density open cosmology (Cl m = 0.3, 
center) and a fiat, matter-dominated model (Cl m = 1, right). 



First of all, we may note that the expected variation with 
cosmic time of the total galaxy luminosity due to evolution in 
its stellar component does not affect the metric definition of 
angular diameters unless this luminosity change depends on 
radius (see paper III of this series (marinoni et al. 2007) for 
a detailed analysis of this issue). Moreover, we can check each 
galaxy spectrum or image for peculiarities indicating possible 
evolution or instability of the standard rod which may be in- 
duced by environmental effects, interactions or excess of star 
formation. 



Any possible size evolution of the standard rod needs to be 
taken into account next. Interestingly, it has been shown by dif- 
ferent authors that large discs in high redshift samples evolve 
much less in size than in luminosity in the redshift range < 
Z < 1. Recent studies show that the amount of evolution to z~l 
appears to be somewhat smaller than expected: disc sizes at 
z ~ 1 are typically only slightly smaller than sizes measured lo- 
cally dTakamival 119991 iFaber et all Exill iNelson et all 12002 ' 



Totani et alll2002h~ This is also theor etically predicted by sim- 



ulations; Bois sier & Prantzosl (1200 II) . for example, show that 
large discs (i.e. fast rotators) should have basically completed 
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Fig. 7. 1,2, and 3<x confidence level contours in the [Qq, w] plane computed by applying the angular-diameter test to data unaf- 
fected (upper panel) and affected by diameter evolution (lower panel). We consider a linear model for disc evolution normalized 
by assuming that discs were smaller by 30% at z — 1.5, and a nominal relative error in the standard rod measures of 5% per 
redshift bin (see discussion in §6.1). The effects of disc evolution onto cosmological parameter estimation are compared to 
the evolution-free case for three different fiducial cosmologies: a low matter density open cosmology (Q,„ = 0.3, left) a fiat, 
A-dominated cosmology (Qa = 0.7, center), and an Einstein-de Sitter model (Q,„ = 1, right). 



their evolution already by z ~ 1 and undergo very little increase 
in size afterwards. Infall models (e.g., Chiappini, Matteucci, & 
Gratton 1997, Ferguson & Clarke 2001, Bouwens & Silk2002) 
also predict a mild disk size evolution. Disk sizes at z ~ 1 in 
these models are typically only 20% smaller than at z = 0. 

6.1. Analysis of the biases introduced by evolution 

Even if literature evidences are encouraging, we have to be 
aware that even a small amount of evolution may introduce 
artificial features and bias the reliability of the cosmological 
inferences. 

In this section we directly address this issue by consider- 
ing different evolutionary patterns for the standard rods, and 
by analyzing to what level the simulated true cosmological 
model may still be correctly inferred using evolved data. In 
other words, we investigate how different disc evolutionary his- 



tories affect the determination of cosmological parameters by 
answering to the following three questions: 

a) is there a feature that may be used to discriminate the 
presence of evolution in the data? 

b) which cosmological parameter is more sensitive to the 
eventual presence of disc evolution?, and in particular what are 
the effects of evolution on the value of the inferred dark energy 
density parameter Og? 

c) is there a particular evolutionary scenario for which the 
inferred values of Q.q and w are minimally biased? 

For the purposes of this analysis, we consider for the an- 
gular diameter-redshift test the baseline (0.1 < z < 1.4) di- 
vided in bins of width dz = 0.1 and assume a relative scat- 
ter in the mean size per bin of 5%. This scatter nearly corre- 
sponds to that expected for a sample of 1300 rotators (with 
< z < 1.4 and dz = 0.1) whose diameters are individually 
(and locally) calibrated with a 40% precision. We then select 
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a given fiducial cosmology (input cosmology), apply an arbi- 
trary evolution to the standard rods, and then fit the evolved 
data with the unevolved theoretical prediction given in eq.[3]in 
order to obtain the best fitting (biased) output cosmology and 
the associated confidence levels contours. We decide that the 
best fitting cosmological model offers a good fit to the evolved 
data if the probability of a worse x 1 1S smaller that 5% (i.e. 

P(* 2 >x 2 obs )<o.05). 

We adopt three different parameterizations to describe an 
eventual redshift evolution of the velocity selected sample of 
galaxy discs D v : a late-epoch evolutionary scenario (AD/D = 
(D v (z) - D v ,(0))/D v (0) = -|(5i| Vz) where most of the evolution 
is expected to happen at low redshifts and levels off at greater 
distances (6\ is the relative disc evolution at z=l), a linear evo- 
lutionary scenario (AD/D = -|5i|z) without any preferred scale 
where major evolutionary phenomena take place (i.e. the gra- 
dient of the evolution is nearly constant), and an early-epoch 
evolution scenario (AD/D - -\6\\z 2 ) where most of the evolu- 
tion is expected to happen at high redshift. 

We note that for modest disc evolution, the linear parame- 
terization satisfactorily describes the whole class of evolution- 
ary models whose series expansion may be linearly represented 
(for example, the hyperbolic model (D v (z) = D v (0)/(1 + \6\\z))). 
For z « 1 it also represents fairly well the exponential model 
(D v (z) = D V (0)(1 - z) s ). Moreover, the linear model is the fa- 
vored scenario for disc size evolution at least at low redshift 



(< 1.5) as predicted by sim ulations (e.g. iMo. Mao & W hite. 
19981; iBouwens & Silkll2002b . 

First, let's assume that w = — 1 and that the dark energy be- 
haves like Einstein's cosmological constant. In Fig.|6]we con- 
sider three different input fiducial cosmological models (a flat 
A-dominated universe (Q A = 0.7), an open model (Q m = 0.3) 
and an Einstein-de Sitter universe (Q m = 1)) and show the char- 
acteristic pattern traced by the best fitting output values (Q,„, 
Qa) inferred by applying the angular diameter test to data af- 
fected by evolution. A common feature of all the various evo- 
lutionary schemes considered is that the value of Q,„ is system- 
atically underestimated with respect to its true input value: the 
stronger the evolution in diameter and the smaller Q m will be, 
irrespectively of the particular disc evolutionary model consid- 
ered. Since many independent observations consistently indi- 
cate the existence of a lower bound for the value of the normal- 
ized matter density (Q m > 0.2), we can thus use this parameter 
as a sensitive indicator of evolution. 

Once the presence of evolution is recognized, the remain- 
ing problem is to determine the level of bias introduced in the 
dark energy determination. If the gradient of the disk evolution 
function increases with redshift (quadratic evolution), then the 
estimates of Qa are systematically biased low. The contrary 
happens if the evolutionary gradient decreases as a function of 
look-back time (square root model). If the disc evolution rate 
is constant (linear model), then even if discs are smaller by a 
factor as large as 40% at z — 1 .5 the estimate of the dark energy 
parameter is only minimally biased. The net effect of a linear 
evolution is to approximately shift the best fitting Qa value in 
a direction parallel to the Q m axis in the [Q m , Qa] plane. 

More generally, by linearly evolving disc sizes so that they 
are up to 40% smaller at z=1.5, and by simulating the apparent 



angle observed in any arbitrary cosmological model with mat- 
ter and energy density parameters in the range < Q m < 1 and 
< Qa ^ 1 (w — -1) we conclude that the maximum deviation 
of the inferred biased value of Qa from its true input value, is 
limited to be \max(6Q.\)\ < 0.2, whatever the true input value of 
the energy density parameter is. In other terms, in the particular 
case of a linear and substantial (< 40% at z=1.5) evolution of 
galaxy discs, the central value of the dark energy parameter is 
minimally biased for any fiducial input model with < Q,„ < 1 
andO < Q A < 1). 

We have shown that the presence of evolution is unambigu- 
ously indicated by the "unphysical" best fitting value of the 
parameter Q,„. We now investigate the amplitude of the biases 
induced by disc evolution in the [Qg, w] plane. We assume for 
this purpose that w is free to vary in the range -1 < w < — 1/3, 
which means assuming that the late epoch acceleration of the 
universe might be explained in terms of a slow rolling scalar 
field. 

We first consider a situation where the disc size evolution 
is modest, and could be represented by any of the three mod- 
els considered. Whatever the mild evolution model considered 
(less then 15% evolution from z=1.5) and assuming a scatter 
in the angular diameter-redshift diagram of 5% in each red- 
shift bin, we find that the input values of Qg and w are con- 
tained within the l<x biased confidence contour derived from 
the evolved data. 

Fig.|7]shows the 1, 2 and 3<x "biased" confidence contours 
obtained by fitting with eq.[3]a simulated angular-diameter red- 
shift diagram in which discs have been linearly evolved. We 
show that, if disc evolution depends linearly on redshift and 
causes galaxy dimensions to be up to 30% smaller at z — 1.5, 
the true fiducial input values of Q.q and w are still within l<x 
of the biased confidence contours inferred in presence of a 
standard rod evolution (and a scatter in the angular diameter- 
redshift relation as low as 5% in each redshift bin). We have 
tested that these conclusions hold true for every fiducial in- 
put cosmology with parameters in the range < Q m < 1, 
< Q G < 1 and -1 < w < -1/3. 

Thus, if, disc evolution is linear (as predicted by theoretical 
models) and substantial (up to ~ 30% at z = 1.5), or arbi- 
trary and mild (up to ~ 15% at z = 1.5), then in both cases 
the angular diameter test reduces from a test of the whole set 
of cosmological parameters, to a direct and fully geometrical 
test of the parameters subset (Qg,w). For example, in a mini- 
mal approach, the angular diameter test could be used to test 
in a purely geometrical way the null hypothesis that "a dark 
energy component with a constant equation of state parame- 
ter w is dominating the present day dynamics of the universe" . 
Moreover, as Fig.|7]shows, an universe dominated by dark en- 
ergy may be satisfactorily discriminated from a matter domi- 
nated universe (Qg = 0). 

In the evolutionary pictures considered, galaxy discs are 
supposed to decrease monotonically in size in the past. Since 
the sensitivity of the test to changes in the linear diameter D is 
described by a growing monotonic function in the redshift in- 
terval [0,1.4], one may hope to test cosmology in a way which 
is less dependent on systematic biases by limiting the sample 
at z < 1 • However, by doing this we would halve the number 
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of standard rods available for the analysis (~ 600 with respect 
to the original ~ 1300). The test efficiency in constraining cos- 
mological parameters would consequently be significantly de- 
graded. 

6.2. The Hubble-Diagram using galaxies 

The same velocity criterion that allows the selection of standard 
rods also allows the selection of a sample of standard candles. 
Using the same set of tracers, the Tully & Fisher relationships 
connecting galactic rotation velocities to luminosities and sizes 
offer the interesting possibility of implementing two different 
cosmological tests, the angular-diameter test and the Hubble 
diagram. 

Thus, we have analyzed how the effects of luminosity evo- 
lution may bias the estimation of cosmological parameters in 
the case of the Hubble diagram. Assuming that the absolute lu- 
minosity L of the standard candle increases as a function of red- 
shift according to the square root, linear and quadratic scenar- 
ios, then the value of O b , is systematically overestimated. The 
value of £2„, is biased in an opposite sense with respect to the 
angular diameter test (see Fig. |6j. Therefore, it is less straight- 
forward to discriminate the eventual presence of evolution in 
the standard candle on the basis of the simple requirement that 
any "physical" matter density parameter is characterized by a 
positive lower bound. The different Q m shifts (with respect to 
the fiducial value) observed when the evolved data are fitted 
using the Hubble diagram or the angualr diameter test are due 
to the fact that, given the observed magnitudes and apparent 
angles, an increase with redshift in the standard candle abso- 
lute luminosity causes the best fitting distances to be biased 
towards higher values, while a decrement in the physical size 
implies that real cosmological distances are underestimated. 

Moreover, even considering a linear evolutionary picture 
for the absolute luminosity as well as a modest change in the 
standard candle luminosity, i.e. Am = M(z) - M(0) = -0.5 at 
z=1.5, the input fiducial cosmology falls outside the 3cr con- 
fidence contour obtained by applying the magnitude-redshift 
test to the sample of evolved standard candles. Note that error 
contours are derived by assuming a scatter of cr M = 0.05 per 
redshift bin in the Hubble diagram. Since, at variance with the 
size of large discs, galactic luminosity is expected to evolve 
substantially with redshift (wit hin the VVDS su rvey, Am 1 



in the I band for M* galaxies (lllbert et all 12005)) we conclude 



that the direct implementation of the Hubble diagram test as 
a minimial test for the parameter subset [CIq, w] using galaxy 
rotation as the standard candle indicator is more problematic. 
As an additional problem, we note that galaxy luminosity is 
seriously affected by uncertainties in internal absorption cor- 
rections, and that due to K-correction, the implementation of 
the Hubble diagram requires multi band images to properly de- 
scribe the rest frame emission properties of galaxies. 

6.3. The cosmology-evolution diagram 

In this section we want to address the more general question of 
whether it is possible to infer cosmological information know- 



ing a-priori only the upper limit value for disc evolution at 
some reference redshift (for example, the maximum redshift 
surveyed by a given sample of rotators). In other terms, we ex- 
plore the possibility of probing geometrically the cosmological 
parameter space in a way which is independent of the specific 
evolution function with which disc sizes change as a function 
of time. The only external prior is the knowledge of an upper 
limit for the amplitude of disc evolution at some past epoch. 

Given an arbitrary model specified by a set of cosmological 
parameters p, and given the observable 0° bs (z), i.e. the apparent 
angle subtended by a sample of velocity selected galaxies (with 
locally calibrated diameters D v (0)), then 



e B (z,p) = GP bs (z)-er(D v (0),z,p) 



(7) 



is the function which describes the redshift evolution of the 
standard rod, i.e eg = (D v (z) - D v (Q))/dg in the selected cos- 
mology. 

Let's suppose that we know the lower and upper limits 
of the relative (adimensional) standard rod evolution (6(z) = 
d (0) ) at some specific redshift z. 

Assuming this prior, we can solve for the set of cosmo- 
logical parameters (i.e. points p of the cosmological parameter 
space) which satisfy the condition 



6,(1) < 



We(p,z)\ 



< S u (l). 



(8) 



ff h (D(0),z,p) 

This inequality establishes a mapping between cosmology 
and the amount of disc evolution at a given redshift which 
is compatible with the observed data. By solving it, one can 
construct a self-consistent cosmology-evolution plane where 
to any given range of disc size evolution at I corresponds in 
a unique way a specific region of the cosmological parame- 
ter space. Vice versa, for any given cosmology one can ex- 
tract information about disc evolution. Clearly, the scatter in 
the angular-diameter diagram directly translates in the uncer- 
tainties associated to the evolution boundaries in the cosmo- 
logical space. 

We note that the boundaries of the region of the cosmologi- 
cal parameter space which is compatible with the assumed prior 
on the evolution of diameters at z can be equivalently expressed 
in term of the maximum absolute evolution in luminosity (see 
appendix A). This because, as stated in the previous section, 
velocity selected objects have the unique property of being at 
the same time standards of reference both in size and luminos- 
ity. 

With this approach, one may by-pass the lack of knowledge 
about of the particular evolutionary track of disc scalelengths 
and luminosities and try to extract information about cosmol- 
ogy/evolution by giving as a prior only the fractional evolution 
in diameters or the absolute evolution in magnitude expected 
at a given redshift. The essence of the method is as follows: 
instead of directly putting constraints in the cosmological pa- 
rameters space by mean of cosmological probes, we study how 
bounded regions in the evolutionary plane (^, AM) map onto 
the cosmological parameter space. 

In Fig [8] we show the cosmology-evolution diagram de- 
rived by solving eq. [8] for different ranges of 6(z). The refer- 
ence model is the concordanace model (D. m = 0.3, SIq = 0.7, 
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w — - 1 ) and the reference evolution at z — 1 .5 is assumed to be 
6(z) = 0.25 for discs, and M v (f) - M v (0)=-1.5 for luminosities. 
The cosmology-evolution diagram represents the unique corre- 
spondance between all the possible cosmological models and 
the amount of evolution in size and luminosity which is com- 
patible with the observed data at the given reference redshift. 
Using independent informations about the range of evolution 
expected in the structural parameters of galaxies, for exemple 
from simulations or theoretical models, one may constrain the 
value of cosmological parameters Vice versa, if the cosmolog- 
ical model is known, then one may directly determine the evo- 
lution in magnitude and size of the velocity selected sample of 
rotators. 

With this approach, the possibility of discriminating be- 
tween different cosmologies depends on the amount of evolu- 
tion affecting the standard rods at redshift z. Two different sets 
(pi and pi) of cosmological parameters may be discriminated 
at z = z if the relative disc evolution at z is known to a precision 
better than 



5(1) < 



\r(z,p2)-r(z,pi)\ 
r(z,pi) 



(9) 



For example, this kind of analysis shows that an Einstein- 
de Sitter universe may be unambigously discriminated from a 
critical universe with parameter Q. m = 0.3, Qg = 0.7, w - — l.if 
the relative disc evolution at z = 1 .5 is known with a precision 
better than 28%. In a similar way, an open (Q m = 0.3, Qg = 
0) universe may be discriminated from an Einstein de-Sitter 
universe if S(z) < 17%. 



7. Conclusions 

The scaling of the apparent angular diameter of galaxies with 
redshift 6(z) is a powerful discriminator of cosmological mod- 
els. The goal of this paper is to explore the potentiality of a new 
observational implementation of the classical angular-diameter 
test and to study its performances and limitations. 

We propose to use the velocity-diameter relationship, cal- 
ibrated using the [OII]/13727A line- widths, as a tool to select 
standard rods and probe world models. As for other purely ge- 
ometrical test of cosmology, a fair sampling of the galaxy pop- 
ulation is not required. It is however imperative to have high 
quality measurements of the structural parameters of high red- 
shift galaxies (disc sizes and rotational velocity). Surveys with 
HST imaging and high enough spectral resolution will thus 
provide the fundamental ingredients for the practical realiza- 
tion of the recipe we have presented. 

In order to avoid any luminosity dependent selection ef- 
fect (such as for exemple Malmquist bias) it is necessary to 
apply the proposed test to high velocity rotators. We show that 
nearly 1300 standard rods with rotational velocity in the bin 
V ~ 200 + 20 km s _1 ) are expected in a field of size 1 deg 2 over 
the redshift range < z < 1 .4. Interestingly this large sample 
can be quickly assembled by the currently underway zCOS- 
MOS deep redshift survey, which uses the VIMOS multi-object 
spectrograph at the VLT to target galaxies photometrically se- 
lected using high-resolution ACS images. 
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Fig. 8. Cosmology-evolution diagram for simulated data which 
are affected by evolution. Apparent angles and luminosities 
of the velocity selected sample of rotators are simulated in a 
ACDM cosmology. Standard rods and candles have been arti- 
ficially evolved so that at I = 1.5 discs are 26% smaller and 
luminosities 1 .4 mag brigther. The cosmological plane is parti- 
tioned with different boundaries obtained by solving equation 
[8]for different values of 6(z = 1 .5), i.e. of the external prior rep- 
resenting the guessed upper limit of the relative disc evolution 
at the maximum redshift covered by data. The external prior 
is also expressed in term of absolute luminosity evolution (see 
discussion in Appendix A). 

Even allowing a scatter of 40% for the 
[OII]/13727A linewidth-diameter relationship for disc 
galaxies, we show that the angular-diameter diagram con- 
structed using this sample is affected by a scatter of only ~ 5% 
per redshift bin of amplitude dz — 0.1. This scatter translates 
into a 20% precision in the "geometric" measurement of the 
dark energy constant equation of state parameter w, through a 
test performed without priors in the [Q. m , Qq] space. 

Current theoretical models suggest that large discs (i.e. fast 
rotators) evolve weakly with cosmic time from z = 1 .5 down to 
the present epoch. Anyway, we have explored how an eventual 
evolution of the velocity-selected standard rods might affect the 
implementation of the test. We have shown that any possible 
evolution in the standard rods may be unanbiguously revealed 
by the fact that even a small decrement with redshift of the disc 
sizes shifts the inferred value of the matter density parameter 
into "a-priori excluded" regions (Q,„ < 0.2). 

We have shown that a linear (as expected on the basis of 
various theoretical models) and substantial (up to 40% over the 
range < z < 1 .5) disc evolution minimally biases the inferred 
value of a dark energy component that behaves like Einstein's 
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Fig. 9. As in fig. [8]but for a different cosmology (an Einstein- 
de Sitter universe). 



cosmological constant A. Moreover we have shown that as- 
suming that discs evolve in a linear-like way as a function of 
redshift, and that their sizes were not more than 30% smaller 
at z — 1 5 with respect to their present epoch dimension, then 
the angular diameter test can be used to place interesting con- 
straints in the [Qg, w] plane. In particular, assuming a scatter 
of 5% per redshift bin in the angular diameter-redshift diagram 
(nearly corresponding to the scatter expected for a sample of 
1300 rotators with < z < 1.4, dz = 0.1, whose diameter is lo- 
cally calibrated with a 40% precision), we have shown that the 
input fiducial [Q.q, w] point is still within the lo - error contours 
obtained by applying the angular diameter test to the evolved 
data. 

Finally, we have outlined the strategy to derive a 
cosmology-evolution diagram with which it is possible to es- 
tablish an interesting mapping between different cosmological 
models and the amount of galaxy disc/luminosity evolution ex- 
pected at a given redshift. The construction of this diagram 
does not require an a-priori knowledge of the particular func- 
tional form of the galaxy size/luminosity evolution. By reading 
this diagram, one can infer cosmological information once a 
theoretical prior on disc or luminosity evolution at a given red- 
shift is assumed. In particular if the amplitude of the relative 
disc evolution at I = 1.5 is known to better than ~ 30%, then 
an Einstein-de Sitter universe (£2 m =1) may be geometrically 
discriminated from a flat, vacuum dominated one (Q m = 0.3, 
Qq = 0.7). Viceversa, one can use the cosmology-evolution 
diagram to place constraints on the amplitude of the galaxy 
disc/luminosity evolution, once a preferred cosmology is cho- 
sen. 

In conclusion, given the simple ingredients entering the 
proposed implementation strategy, nothing, besides evolution 
of discs, could in principle bias the test. Even so, evolution can 



be easily diagnosed and, under some general conditions, it can 
be shown that it does not compromise the possibility of detect- 
ing the presence of dark energy and constraining the value of 
its equation of state. 



In the following paper s of this series (ISaintonge et al 



2007; M arinoni et aUl2007h . we implement the proposed strat- 



egy to a preliminary sample of velocity-selected high redshift 
rotators. 
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8. Appendix A 

Given a spectroscopically selected sample of objects with con- 
stant rotational velocity we can derive the observed magnitude 
m° of a standard candle of absolute mag nitude M v ( 0) loc ated at 
redshift z, by using the standard relation Sandage d 19721) 



m° = m' h (M v (Q),z,p) + e M (z) + K(z) (10) 
where 

m th = M v (0) + 5 logd L (z,p) + 25 

and were di is the luminosity distance, K(z) is the K correction 
term and 6m(z) is the a-priori unknown evolution in luminosity 
of our standard candle, i.e. e M (z) = AM v (z) = M v (z) - M v (0) is 
the difference between the absolute magnitude of an object of 
rotational velocity V measured at redshift z and the un-evolved 
local standard value M v (0). 

From the definition of wavelength-specific surface bright- 
ness p. we deduce that the variation as a function of redshift in 
the average intrinsic surface brightness (within a radius R) for 
our set of homologous galaxies is 



A^" , (z)> s = AM„(<«)-51og 



R(z) 
R(0) 



By opportunely choosing the half light radius D v as a metric 
definition for the size of a galaxy we immediately obtain 

A(/-( Z ))^^, P )-51og( ^ ( ^ p) + l) (11) 

The intrinsic surface brightness evolution is not an observ- 
able, but in a FRW metric this quantity is related to the surface 
brightness change observed in a waveband A/l by the relation 



A(p°(z))d = Hp"\z)) D + 2.5 log(l + z) 4 + K(z) 



(12) 



Thus, once we measure the redshift evolution of A(p°(z))d 
for the sample of rotators, the absolute evolution in luminosity 
corresponding to a given relative evolution in diameters can be 
directly inferred using eq. Qj] 



Marinoni et al.: Geometrical tests of cosmological models. I, 



15 



References 

Aldering, G. et al. 2002, in SPIE Proceedings Vol 4835, 

|astro-ph/0209550| 
Berlind, A. A., & Weinberg, D. H. 2002, ApJ, 575, 587 
Boissier, S., Prantzos, N. 2001, MNRAS, 325, 321 
Bottinelli, L., Gouguenheim, L., Paturel, G, & de Vaucouleurs, 

G. 1980, ApJ, 242, L153 
Bouwens, R., & Silk, J. 2000, ApJ, 568, 522 
Bruzual, A. G, & Spinrad, H. 1978, ApJ, 220, 1 
Buchalte, A., Helfand, D. J., Becker, R. H., & White, R. L. 

1998, ApJ, 494, 503 
Carroll, S. M., Press, W. H., & Turner, E. L. 1992, ARA&A, 

30, 499 

Chiappini, C, Matteucci, E, & Gratton, R. 1997, ApJ, 477, 765 
Colless, M., Glazebrook, K., Mallen-Ornelas, G, & 

Broadhurst, T. 1998, ApJ, 500, L75 
Cooray, A., Hu, W., Huterer, D., & Joffre, M. 2001, ApJ, 557L, 

7 

Daly, R. A. 2004, ApJ, 612, 652 

Daly, R. A., & Djorgovski, S. G. 2004, ApJ, 612, 652 

Davis, M, Newman, J. A., Faber, S. M, & Phillips, A. C. 

2000 in Deep Fields, Proceedings of the ESO/ECF/STScI 

Workshop, eds. Cristiani, S., Renzini, A., Williams, R. E., 

Springer, p. 241 
de Bernardis, P. et al. 2002, ApJ, 564, 559 
de Vaucouleurs, G, de Vaucouleurs, G, & Corwin, H. G. 1976, 

Second Reference Cataloguwe of Bright Galaxies( Austin: 

University of Texas) RC2 
Djorgovski, S., & Spinrad, H. 1981, ApJ, 251, 417 
Faber, S. M., Phillips, A. C, Simard, L., Vogt, N. P., & 

Somerville, R. S., 2001, in Galaxy Disks and Disk Galaxies, 

ASP Conference Series, Vol. 230, p. 517 
Ferguson, A. M. N., & Clarke, C. J. 2001, MNRAS, 325, 781 
Gonzalez, A. H. 2002, ApJ, 567, 144 

Gurvits, L. I., Kellermann, K. I., & Frey, S. 1999, A&A, 342, 
378 

, Haiman, Z., Mohr, J. J., & Holder, G. P. 2001, ApJ, 553, 545 
Halverson, N. W., et al. 2001, ApJ, 568, 38 
Hickson, P., 1977a, ApJ, 217, 16 
Hogg, D. W. 1999, |astro-ph/9905116| 

Hoyle, F, 1959 in Proc. IAU Symp. 9 and URSI Symp. 1, Paris 
Symp on Radio Astronomy, ed. R. N. Brachewell (Stanford: 
Stanford Univ. Press), 529 

Hubble, E., & Tolman, R. C, Astrophys. J., 1935, 82, 302 

Huterer, D., & Turner, M. S. 2001, Phys.Rev. D64 

Ilbert, O., et al. 2005, A&A, 439, 863 

Kapahi, V. K., 1975, MNRAS, 172, 513 

Kellermann, K. I. 1993, nature 361, 134 

Kobulnicky, H. A., & Gebhardt, K., ApJ, 1999, 119, 1608 

Lahav, O. 2002, MNRAS, 333, 961 

Lee, A. T. et al. 2001, ApJ, 561, LI 

Le Fevre, 0. et al. 2005, A&A, 439, 845 

Lilly, S. J., Schade, D., Ellis, R., Le Fevre, O., Brinchmann, J., 
Tresse, L., Abraham, R., Hammer, F, Crampton, D., Colless, 
M., Glazebrook, K., Mallen-Ornelas, G, & Broadhurst, T. 
1998, ApJ, 500, L75 

LiUy, S. J., et al. 2006, ApJ, in press ( |astro-ph/061229T] l 



Lima, J. A. S., & Alcaniz, J. S. 2002, A&A, 566, 15L 
Marinoni, C, Monaco, P., Giuricin, G, & Costantini, B. 1998, 

ApJ, 505, 484 
Marinoni, C. & Hudson, M. 2002, ApJ, 569, 101 
Marinoni, C, et al. 2006, A&A, 442, 801 
Marinoni, C. et al. 2007, A&A in press, astro-ph/07 10.0761 
Miley, G. K., 1971, MNRAS, 152, 477 

Mo, H. J., Mao, S., & White, S. D. M. 1998, MNRAS, 295 319 
Narayan, R., & White, S. D. M. 1988, MNRAS, 231, 97 
Newman, J, A., Marinoni, C, Coil, A. L., & Davis, M. 2002, 

PASP, 1 14, 29 
Newman, J. A., & Davis, M. 2000, ApJ, 534, LI 1 
Nelson, A, E., Simard, L., Zaritsky, D., Dalcanton, J. J., & 

Gonzalez, A. H., 2002, ApJ, 567, 144 
Pen, U-L. 1997, New A., 2, 309 
Perlmutter, S. et al. 1999, ApJ, 517, 565 
Petrosian, V. 1976, ApJ, 209, LI 
Riess, A. G. et al. 1999, ApJ, 560, 49 

Saintonge, A., Masters, K.L., Marinoni, C, Giovanelli, R., & 

Haynes, M.P 2007, A&A submitted, astro-ph/07 10.0760 
Sandage, A., ApJ, 1972, 173, 485 

Sandage, A., 1998, Ann. Rev. Astron. Astrophys. 26, 561 
Sandage, A., 1995, in Saas-Fee Advanced Course 23, The Deep 

Universe: Practical Cosmology; ed. B. Binggeli & R. Buser 

(New York: Springer) 
Spergel, D. N., et al. 2006, |astro-ph/0603449| 
Springob, C. M., Masters, K. L., Haynes, M. P., Giovanelli, R„ 

& Marinoni, C. 2007, ApJS in press, |astro-ph/0705064| 7 
Takamiya, M., 1999 ApJS, 122, 109 
Tegmark, M. 2006, Phys. Rev. D 74, 123507 
Tolman, R. C, 1930, Proc. Natl. Acad. Sci. USA: 16, 511 
Totani, T, Yoshii, Y., Maihara, T, Iwamuro, F, & Motohara, 

K. 2001, ApJ, 559, 592 
Tully, R. B., & Fisher, J. R. 1977, A&A, 54, 661 
Turner, M. S., & White, M. 1997, Phys. Rev. D, 56, 4439 
van den Bosch, F. C, et al. 2006, MNRAS submitted 

( |astro-ph 0610 686) 
Wilkinson, P.N., Browne, I. W. A., Alcock, D. et al. 1998 in 

Observational Cosmology with new radio surveys, Kluwer, 

Acad. Publishers, Dordrecht, p. 221 



